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1.  Introduction  1*  1  i 

Satellites  of  the  NNSS  navigation  system  (Navy  Naviga¬ 
tion  Satellite  System),  or  TRANSIT  satellites,  which  in  1967 
were  made  available  for  commercial  flight  navigation,  have 
found  widespread  application  in  geodetic  work.  At  present, 
many  methods  utilize  NNSS  satellites  to  determine  the  coord¬ 
inates  of  the  observer  on  the  basis  of  Doppler  measurements 
and  the  known  position  of  satellites  [1,  5,  6], 

The  position  of  the  SN  (navigational  satellites)  in 
the  system  considered  Is  calculated  on  the  basis  of  data 
transmitted  in  2-mlnute  time  Intervals  from  the  satellite 
memory  block  located  on  the  ship.  During  this  time  the  SN 
moves  about  900  km  in  space,  which  allows  determination  of 
the  coordinates  of  the  observer  from  Doppler  measurements. 

The  accuracy  of  determination  of  coordinates  depends,  among 
other  things,  on  geometric  conditions  of  the  SN  trajectory 
above  the  horizon  of  the  observer.  The  Doppler  measurements 
made  during  the  most  favorable  flights  of  the  SN  above  the 
horizon  are  used  In  calculations.  For  an  effective  utiliza¬ 
tion  of  measurements,  and  for  shortening  the  computation  time, 
it  is  necessary  to  have  information  about  the  SN  trajectory. 


In  other  words,  it  is  necessary  to  know  their  ephemerides. 

The  knowledge  of  ephemeris  is  indispensable  for  turning 
on  the  measurement  instruments  for  receiving  signals  from 
the  SN.  This  knowledge  is  also  necessary  for  receivers  which 
automatically  receive  signals  from  the  satellite  as  soon 
as  it  appears  over  the  horizon.  If  there  are  overflights  of 
two  SN  at  approximately  the  same  time,  the  receiver  will 
accept  signals  from  the  satellite  which  appeared  first  over 
the  horizon,  and  will  ignore  signals  of  the  second  satellite. 
Only  an  operator  can  switch  the  SN  signal  receiver  to  a  signal 
with  more  favorable  geometric  conditions.  On  the  basis  of 
data  transmitted  from  the  SN,  one  can  calculate  the  ephemeris 
for  any  point  on  the  Earth  many  weeks  ahead. 

This  paper  gives  practical  formulas,  tables  and  graphs 
which  enable  fast  determination  of  flight  conditions  of  naviga¬ 
tional  satellites  of  the  NNSS  system,  particularly  for  the 
terrain  of  Poland.  This  is  preceded  by  an  example  of  decipher¬ 
ing  Information  sent  by  a  SN,  and  calculation  of  its  position 
In  space  on  the  basis  of  transmitted  data. 

2.  Decloherinc  of  NNSS  Message 

NNSS  navigation  satellites  have  memory  blocks.  The  antici 
pated  elements  of  the  Kepler  orbit,  the  so-called  orbital 
parameters,  and  information  about  time  is  fed  into  these  memory 
blocks  covering  a  period  of  about  12  hours.  These  data  are 
necessary  to  find  the  position  of  the  SN  in  the  geocentric, 
system  and  are  transmitted  in  two-minute  intervals  on  both 
carrier  frequencies  (150  MHz  and  *J00  MHz)  without  disturbing 
the  Doppler  measurements.  Digital  signals  sent  in  2 -minute 
intervals  contain  6103  binary  bits,  made  up  of  156  words  of  39 
bits  each  and  the  last  157-th  word  of  19  bits.  These  words 
form  25  groups,  each  having  6  words  of  39  bits  (25  groups). 


The  last  group  contains  7  words.  One  group  of  6  words  is 
transmitted  during  4S600908.  This  interval  or  its  multiple 
is  utilized  in  the  short  Doppler  method.  The  words  with 
numbers  1,  2,  3  contain  information  about  time,  and  between 
the  transmission  of  words  2  and  3  there  is  the  so-called  time 
marker  given  with  the  highest  possible  accuracy,  of  the  order 
of  a  few  microseconds.  The  time  marker  refers  to  each  even  TU 
minute,  and  is  the  initial  time  for  the  Doppler  calculations. 

Of  the  remaining  154  words,  only  one  group  of  25  words  (No. 

8,  14,  20,  ...152)  is  utilized  in  navigation  and  geodetic 
practice.  The  information  contained  in  the  remaining  words 
has  not  been  disclosed;  it  is  of  a  military  nature  and  is  not 
deciphered  by  civilian  receivers  [1,2]. 

The  words  of  interest  to  us  (with  numbers  8,  14,  ...152) 
are  nine-digit  numbers  in  the  decimal  system.  The  first  8 
words  (with  numbers  8,  14,  20,  26,  32,  38,  44,  50)  contain 
information  about  the  so-called  orbital  parameters.  The  content 
of  one  word  refers  to  one  TU  expressed  in  even  minutes.  During 
each  transmission  we  are  getting  ephemeris  parameters  for  a 
14-minute  interval  of  time  (every  2  minutes),  which  enables 
their  interpolation  for  the  desired  intermediate  times.  This 
is  done  for  instance,  in  short  Doppler  programs.  Each  ephemeris 
word  gives  the  number  of  2-minute  intervals  of  time  which  have 
elapsed  since  the  last  full  or  half  hour  of  TU  time.  Once  the 
TU  time  is  known  with  an  accuracy  of  several  minutes  or  even 
more  than  ten  minutes  from  the  above  mentioned  data,  we  assign 
a  time  to  the  beginning  of  a  given  transmission,  and  times 
to  particular  words  of  the  ephemeris.  Moreover,  in  each  ephem¬ 
eris  word  there  is  an  actual  correction  to  the  eccentric  anomaly 
AE,  semimajor  axis  Act,  and  partial  information  about  the 
component  nonral  to  the  Keplerian  orbit  T. 

The  deciphering  of  ephemeris  words  is  rather  complex.  The 
first  digit  of  the  ephemeris  word  contains  information  about 
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the  first  time  information  digit  nfcJ  and  information  about 
the  signs  of  corrections  AE  and  Aot.  The  deciphering  of  the  /212 
first  ephemeris  word  digit  n  is  given  in  Table  1  (the  first 
four  columns).  Now,  from  the  second  ephemeris  word  digit 
we  obtain  the  number  of  two-minute  intervals  which  have 
elapsed  since  the  last  full  hour  (half  hour)  TU. 


Table  1. 
("Znak  *  sign) 


The  next  3  digits  give  the  correction  AE  in  units  of  0°0001. 
The  digits  in  positions  6,  7  and  8  of  the  ephemeris  word  give 
the  value  of  the  absolute  correction  Aot  in  units  of  0.01  km. 
The  last  (nineth)  digit  of  the  ephemeris  word  contains  informa 
tion  about  the  component  Z.  In  order  to  decipher  It  properly, 
one  has  to  know  accurately  the  time  to  which  the  given  ephe¬ 
meris  word  refers.  The  first  digit  of  Z  is  determined  on  the 
basis  of  the  value  of  the  9th  digit  of  the  ephemeris  word  n^ 
in  relation  to  TU  expressed  In  minutes  and  divided  by  4. 

Its  deciphering  is  given  In  the  last  two  columns  of  Table  1. 
The  second  digit  of  component  1  is  left  unchanged,  as  the 
ninth  digit  of  the  next  ephemeris  word.  It  is  written  next 

to  the  value  Y  taken  from  Table  1.  Hence,  the  value  of 
P 

component  1  is  given  every  4  minutes.  This  value  consists  of 
two  digits  and  a  sign,  and  it  expresses  Z  in  units  of  0.01  km. 

After  each  consecutive  2-minute  transmission,  the  values 
of  particular  ephemeric  words  move  downwards  by  one  line.  For 
instance,  the  word  number  44  goes  into  the  word  No.  50,  and 
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the  word  No.  50  from  previous  transmission  is  dropped.  The 
location  of  word  No.  8  from  a  previous  transmission  is  now 
taken  over  by  a  new  value.  The  parameters  pertaining  to 
the  given  2-minute  transmission  are  coded  in  word  No.  26  (let 
this  time  be  called  t).  Hence,  the  information  contained  in 
word  No.  8  refers  to  the  moment  t  -  6m,  and  correspondingly 
14  -*■  t  -  4m,  20  -*■  t  -  2m,  32  -*■  +  2m,  38  -*■  t  +  4m,  44  -*•  t  +  6m, 
and  50  -*■  t  +  8m. 

Beginning  with  number  56  the  words  contain  the  actual 
elements  of  the  Keplerian  orbit.  On  the  average,  these  values 
undergo  changes  every  12  hours.  The  meaning  of  particular 
words  is  as  follows: 

the  word  No.  56  contins  value  of  the  time  interval  which 
has  elapsed  from  midnight  GMT  to  the  time  of  the  first 
transit  of  the  satellite  through  perigee  after  the  last 
introduction  of  data  into  the  SN  memory  block  of  (tp); 

-  the  word  No.  62  contains  the  value  of  average  motion  n; 
the  word  No.  68  contains  the  value  of  the  argument  of 
perigee  u  at  the  time  t  ; 

the  word  No.  74  contains  the  absolute  value  of  the  change 
of  argument  of  perigee  jw  | ; 

the  word  No.  80  contains  the  value  of  excentricity  of 
the  orbit  e; 

the  word  No.  86  gives  value  of  the  semimajor  axis  a; 

word  No.  92  gives  value  of  the  longitude  of  the  ascending 

node  ft  at  time  t  : 

P 

word  No.  98  gives  value  of  the  change  of  longitude  of  the 

O 

mode  ft; 

Word  No.  104  gives  value  of  the  cosine  of  the  inclination 
of  the  orbit,  cos  i; 

Word  No.  110  gives  value  of  Greenwich  sidreal  (star)  time 
S  at  the  moment  t  ; 

-  word  No.  116  gives  value  of  the  number  of  the  navigation 
satellite  N; 
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word  No.  122  contains  the  value  t  ,  of  the  time  of 

w 

introducing  data  into  the  SN  memory  block  and  the 
number  of  current  days  counted  from  the  beginning  of 
the  given  year; 

-  word  No.  128  gives  value  of  the  sine  of  the  inclination 
of  the  orbit,  sin  i; 

-  word  No.  134  contains  value  of  Af  -  deviation  from  the 
frequency  400  MHz  transmitted  by  the  SN  transmitter; 
this  value  is  expressed  in  units  of  80  x  10“^; 

words  with  numbers  140,  146  and  152  usually  contain  zeros. 

The  first  decimal  number  of  the  word  No.  56  always  has 
the  value  of  0  or  4,  and  is  deciphered  as  +0,  or  +1  in  the  case 
of  4. 


The  first  digit  of  words  with  numbers  62,  68,  ...134  is 
always  8  or  9;  it  gives  the  sign  of  a  given  value,  the  sign 
plus  (+)  is  denoted  by  8  and  the  sign  minus  (-)  by  9.  The  last, 
i.e.,  the  9th  figure  of  words  62,  68  ...135  is  always  zero  and 
may  be  omitted. 

The  elements  of  the  orbit,  defined  by  fixed  parameters 
contained  in  words  No.  56,  ...»  128,  extrapolate  the  orbit 
of  a  given  SN.  The  elements  of  the  orbit  actually  undergo 
constant  changes  because  of  the  action  of  various  perturbing 
forces,  caused  mainly  by  the  gravitational  field  of  the  Earth 
and  friction  of  upper  atmosphere.  The  effect  of  this  part  of 
the  perturbations  which  can  be  predicted  is  given  in  the 
parameters  of  the  ephemeris. 

Table  2  gives  an  example  of  a  printout  by  a  navigational 
receiver  of  a  message  transmitted  from  SN-18,  and  gives  its 
detailed  deciphering. 

3.  Position  of  SN  in  Space 

On  the  basis  of  data  contained  in  Table  2,  we  shall  calculate 
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the  Greenwich  geocentric  coordinates  of  SN-18.  The  numerical 
example  refers  to  the  time  t  =  13  00  TU  on  17  December  1975. 
Using  formulas  for  Keplerian  motion,  we  calculate  value  of 
the  mean  anomaly  M. 

M  =  n(»  — «j)  =  n(780“  —  446*0124)  w  44°6X5U  (  1  ) 

Taking  into  consideration  this  value  of  M  in  solving  the 
Kepler  equation  Eq  *  M  +  e  sin  EQ  (by  the  method  of  successive 
approximations  or  by  Newton's  method)  we  obtain  the  value  of 
the  excentric  anomaly  Eq  =  H4°9l8ll.  Taking  from  Table  2  the 


Table  2 


Print¬ 

out 


41 414 144  4 
41  011  IMS 


mssuhi 

•IMS  MSS 
as  ms  in* 
•SMS  MSI 
44  4*1 041 4 


Deciphered 

parameters  of  ephemeris 


+raus 


+«.M 

+M4 

7-UW 

TKu 

-•141 

+M» 


"+*•?* 

■+4J4 

+44* 


Orbital  elements  (permanent  part  of  message 


TC 


4  44441144  I,  -  44**4144  -  1 

•  14 71143a  1*4*7HH  ■■  -4. 

4  341411a*  ot  -  147*4114 

•MM144  \U -naamn hm-i 

1  444T4444  .  -  4,447444  ” 

S47441SM  4j.7441.SHas 

5  27442S44  -  ...  174*4114 

4  4*444444  Q  -  4*0044444  Mil -I 
I  00441244  •  >  -  4,044434 

•  19444444  '  S  -  144* 4447 

4  0S4144M  N  m  1414  »■»  SW 
4 14441114  i»-  lMxS"  -ST4»  -  «*MMTV  *81  4s 

iM  17*4.  xn  ms*, 
a  10404444  rts  i  -  1.444444 

4  44171444  —  44,4171  •  14-*  •  444  MBs  —  SS>StT  fcBa 


c*» 

m 


i  ms*.,  u 


Note:  1)  Transmission  from  the  SN  gives  only  the  fractional 
part  of  the  mean  motion.  This  value  has  to  be 
always  preceded  by  the  number  3  denoting  degrees; 

2)  For  SN  TRANSIT  the  rate  of  change  of  the  argument 
of  perigee  (the  abscissa)  always  has  a  negative 
value  (  cii<0,  since  i  =  90°). 

The  deciphering  gives  nominal  accuracies  and  units 
in  which  the  transmitted  values  are  expressed. 

Values  of  a),  (2,  S  refer  to  time  t  . 

P 

correction  AE  *  0°0051  and  adding  it  to  Eq  we  obtain: 

E-J£t  +  J£-  44*92321  'l'. 
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The  value  of  semimajor  axis  is 


7463,78  km  (  3  ) 

In  turn,  we  calculate  the  radius  r  from  the  formula 

t  —  a(i  —  ich£)  =  7424,2024  km  ^  ^ 

In  practice,  we  calculate  first  the  orbital  coordinates  of 
the  SN  (the  system  origin  is  in  the  center  of  Earth,  the  axis 
X  Is  directed  to  perigee,  the  axis  y  lies  in  the  plane  of  the 
orbit  and  is  directed  in  the  direction  of  motion  of  the  SN, 
and  the  axis  ~Z  -  normal  to  the  plane  of  orbit  -  compliments 
the  system  to  make  it  clockwise).  Calculations  are  performed 
with  the  aid  of  formulas 

x  «  a(cotS  —  e)  »  5228,862  km 
y  =■  ay'l—e*  lin E  —  5270,464  kat  (  5  ) 

/213 

2  =  0.08  km  -  as  a  parameter  of  ephemeris  taken  from  Table  2 
Check:  r  =  (  x2  +  y2  +  z2)1//2  =  7424.202  km,  agreeing  with  (4) 

From  the  rectangular  orbital  system  we  can  transform  into 
the  rectangular  equatorial  ( equinocturnal )  or  Greenwich 
system  by  means  of  vectors  F,  H  [4] ,  which  are  functions 
of  the  following  orbital  elements:  ft,  i  ,  ui.  Since  in  practical 
satellite  applications  we  are  interested  in  direct  transforma¬ 
tion  into  the  Greenwich  system,  instead  of  angle  ft  we  substitute 
the  angle  ft  -  S  *  For  a  given  moment  t  we  obtain 

ii  -  S  -  U(tp)  -  >'(!,,)  f-  (U  -k)(l~tp)  *  -  4°  50907 

u>  —  o<( ip)  -  'tui  (i  -  ip)  —  34b  34042  (  6  ) 

o  .  o 

Values  of  the  quantities  ft(t  ),u(t  ),  ft, loo  ,  sin  i  and  cos  i 

W  Q  ^ 

are  obtained  from  Table  2.  S  -  0°250684^76  min-  is  the  value 
of  the  average  rate  of  rotational  motion  of  the  Earth. 

The  components  of  vectors  P,  F  are  calculated  from 
equations 


M2  •  ** SJ? ~£)  ~~ 8in  ®»in<0  —  S)  cos  i  0,9687019 

pUS?  I  JfffrSwm*' °-07MM7 

*  **  «  sin  (O-S)  cos  i  - 

SrAwaW  ,fa  +  -  *■  <®-n  cos  <  - 

Qt  -  cos  o*  sia  I  «  0,9717159 

5*  "  *“  rSf“a2LJln  1  "  —0.0788274 


(7) 


Check:  |  F  j  =  1,  |$|  =  1,  |ff|  =  1,  P- Q  =  F*  R  =  Q.R  =  0. 


On  the  basis  of  (5)  and  (7)  we  obtain  the  desired  SN 

rectangular  Greenwich  coordinates  from  equations 

m  ■  “  6306,1X8  km 

y  -  i*V  +y«,+i«,  -  ~  49647*  km 
.  s-*Pt+yQ,  +  il*i-  388*4»4ta»  0  ' 


Let  r  denote  the  geocentric  radius  of  the  SN,  X  and(f>  geocen¬ 
tric  longitude  and  latitude  -  the  coordinates  of  the  sub¬ 
satellite  point.  When  we  take  into  consideration  the  follow¬ 


ing  formulas 


*  —  rcMfcwl 
y  —  r  cm  9  mb  A 
*  —  rtinp 


(9) 


and  values  from  (8),  we  obtain  y  -  7^24.202  km,  X  =  — il°it988 , 

<p  =  31°56n. 


Knowing  the  position  of  the  SN  in  the  rectangular  geo¬ 
centric  system  and  the  position  of  an  observer  in  the  same 
system,  or  in  a  geodetic  system  on  a  given  reference  ellipsoid, 
one  can  calculate  the  horizontal  coordinates  of  the  satellite. 
We  are  particularly  interested  in  the  horizontal  height 
(altitude)  h .  This  parameter  is  necessary  when  calculating 
corrections  because  of  tropospheric  refraction. 

Predicted  Position  of  NNSS  Satellites 


In  predicting  flight  conditions  of  a  particular  SN  over 
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4. 


a  given  location  on  the  surface  cf  the  Earth,  the  accuracy 
with  which  one  can  extrapolate  changing  elements  of  the  orbit 
over  larger  time  intervals  is  very  important. 

Table  3  gives  transmitted  orbital  elements  of  satellites 
No.s  13,  1^,  12,  19  and  20.  The  elements  of  SN-18  were  given 

in  Table  2.  All  these  data  refer  to  satellites  operating 

in  December  1975. 

Table  4  lists  transmitted  orbital  elements  of  satellites 
No.s  20,  14,  13,  12  and  19  observed  in  January  1976. 

In  Tables  3  and  4,  the  parameters  n,  |w|,  and  ft  are 

expressed  in  units  of  1  degree/minute  angular  velocity.  The 
semimajor  axis  a  is  in  km,  and  the  remaining  parameters  and 
orbital  elements  are  given  in  Tables.  The  second  part  of  Table  4 


Table  3 


/  -SataHu 

Data 

11 

11  XII 

14 

14  Xll 

12 

17  XII 

19 

17  XII 

»  ! 

.7  xn  1 

1 

tp 

S66-9S54 

6  76- SI  11 

235— 52U 

317-1304 

1247-7302  I 

n 

3,3443231 

3,3720142 

5,3815153 

3,3655759 

3.4104,4,  1 

V 

177*1699 

170*5254 

244*8926 

0S°44S6 

1 77*2318  i 

14  1 

0,0019914 
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gives  elements  go  ft  calculated  from  equations 

Lf  C 

®e  =  tu  +  (i»— 

Qe  =  a+liJt  do) 

on  the  basis  of  initial  data  for  w,  ft,  n,  I  to  |  ,  ft  taken  from 

Table  3.  The  calculations  were  made  for  appropriate  times 

t  of  particular  SN  -  data  in  Table  4.  Hence  At  =  t  ,.-t 
P  P3 


*»£•#■*** 


Table  4 
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is  the  time  interval  which  elapsed  between  the  transit  of 
the  SN  over  perigee  given  in  Table  3  and  the  same  moment 

given  in  Table  4.  Table  4  also  gives  the  differences 

(a,  ft 

oj-oig  between  elements^ taken  from  Table  4  and  those  cal¬ 

culated  using  the  values  u>  ,  fi  . 

c  c 


Comparing  the  orbital  elements  of  particular  SN  fron 
Table  3  and  those  after  5-6  weeks  (Table  4),  we  see  that  the 
accurate  calculation  of  predicted  position,  or  ephemeris, 
is  affected  most  significantly  only  by  the  difference  u)-gjc, 
which  in  some  cases  is  smaller  than  1°.  The  remaining 
elements  remain  the  same,  or  change  to  a  small  extent.  The 
sidereal  (star)  time  Sc,  calculated  from  equation 

St  «*  S+$4t 


(11) 

differs  very  little  from  the  corresponding  values  of  S  in 
Table  4,  as  a  result  of  the  simplified  form  of  equation  (11). 


Since  orbits  of  TRANSIT  system  satellites  are  nearly 
circular  and  polar  (e  s  0,  and  i=90°),  in  calculations  of  the 
ephemeris  we  assume  that  e  =  0,  i  =  90°  (the  circular  and  polar 
orbit).  Moreover,  we  assume  that  the  Earth  is  a  sphere  with 
a  mean  radius  R  =  6365  km.  These  simplifying  assumptions 
facilitate  the  calculations  considerably,  and  still  enable  us 
to  obtain  sufficient  accuracy  when  computing  the  ephemerides. 

On  the  basis  of  transmitted  orbital  elements  we  shall 
first  calculate  the  longitude  of  the  ascending  node  X  ,  that 
is  the  geographic  longitude  of  a  point  on  the  equator  at  which 
the  SN  passes  from  the  Southern  to  the  Northern  hemisphere,  /214 
and  also  the  transit  time  through  this  point 

(12) 

The  SN  in  circular  orbit  moves  with  a  uniform  angular  rate  n. 


J 
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and  its  longitude  along  the  orbit  u  at  the  time  t  is  calca 
lated  from  equation 


Let  the  longitude  of  the  node  Aw  be  its  longitude  at  the  first 
transit  of  the  SN  through  the  equator  from  South  to  North  at 
the  time  t  transmitted  in  the  message.  Then  u(t,  )  =  36O0. 
Substituting  this  value  in  (13)  we  get 

360*- •  <9  -  (1*0 

*  _  «  -L. _ LSI 


Let  T  be  the  time  period  required  by  the  SN  to  fly  around 
the  Earth.  It  is  expressed  by  equation 

360*  (15) 

T”  ■— |i| 


It  is  called  the  nodal  period.  After  one  orbit,  the  longitude 
of  the  initial  node  will  change  by  AA^,  and 


therefore  after  k  orbits  we  have 


(16) 


where : 


t° 

w 


(17) 


initial  value  of  the  node,  and  the  transit  time 
of  the  SN  through  this  point,  calculated  from 
equations  (12)  and  (14), 


On  the  basis  of  given  values  of  parameters  A  ,t  we 

w  w 

calculate  spherical  coordinates  of  subsatellite  points  (their 
geocentric  longitude  and  latitude)  at  any  moment  t.  The  long¬ 
itude  is  calculated  from  equation 

*<*>-  a*+(A-ft(«-U 

We  assume  that  the  geographic  longitude  is  calculated 
positive  to  the  East  of  the  Greenwich  meridian  (0°<Ag<  180°), 
and  negative  to  the  West  of  the  Greenwich  meridian 
(»180°<A  <  0°).  Hence,  after  the  SN  passes  through  the  North 
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pole  from  the  Eastern  to  Wesern  hemisphere  we  subtract  180° 
from  the  value  calculated  from  equation  (18).  Conversely, 
we  add  l80°  after  passing  the  pole  from  the  Western  to 
Eastern  hemisphere. 


When  calculating  the  latitude  $  we  have  the  following 
three  cases: 

a)  <J>  =  u,  the  SN  moves  from  the  equator  to  the  North 

Pole 

(19) 

b)  <t>  ~  l80°-u,  the  SN  moves  from  the  North  Pole  to  the 
South  Pole 

for  (20) 

c)  <t>  =  -360°  +  u,  the  SN  moves  from  the  South  Pole  to 
the  equator 

for  <Ua  370*  <  ■  <  3M*  (21) 

where : 

.  (22) 

»m  (»-!•)) (I- tj 


Let  the  coordinates  of  an  observer  (coordinates  of  the 
point  for  which  we  want  to  calculate  the  ephemeris)  be  A0,<J>0. 

First,  we  calculate  the  time  of  apogee  t^.  For  an  SN  in  polar 
orbit,  the  SN  in  at  apogee  is  approximately  at  the  latitude 
of  the  observer.  For  an  observer  in  the  Northern  hemisphere 
two  cases  can  arise,  i.e.,  case  a)for  which  equation  (19) 
applies,  or  case  b)  for  which  equation  (20)  is  relevant  . 


For  apogee,  the  longitude  of  the  satellite  in  orbit  is 


b)  »-  100* -f,' 


(23) 


Knowing  u  we  calculate  the  time  of  apogee  t^  by  substituting 
(23)  into  (22) 
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b)  «{=*.+ 


mu'  -  r. 


(24) 

The  longitude  of  the  subsatellite  point  at  the  apogee  of  the 
SN  is  expressed  by  the  equation 

(25) 


When 

Xk 

^X 0 , 

the 

SN 

moves 

in 

the 

Eastern 

side 

of 

m 

When 

*X0 , 

the 

SN 

passes 

through  zenith 

1 

When 

*  x  0 , 

the 

SN 

moves 

in 

the 

Western 

side 

of 

sky . 

On  the  basis  of  equations  (18)  and  (24)  we  can  calculate 
the  longitude  of  the  ascending  node  for  an  NSSZ  passing  through 
zenith.  The  condition  to  be  fulfilled  is 


4,-  <6 -$><«»-•*> 


(26) 


For  the  cases  a)  and  b)  we  obtain,  respectively 


(27) 

(28) 


The  sign  (-)  refers  to  the  case  in  which  Xo>0°  (the  observer 
is  in  the  Eastern  hemisphere),  and  the  sign  (+)  applies  to 
the  situation  when  Xo<0°  (the  observer  is  in  the  Western 
hemisphere ) . 

In  turn,  we  calculate  the  maximum  interval  of  time  of  the 
SN  stays  over  the  horizon,  which  takes  place  when  a  given 
satellite  passes  through  zenith  (apogee  at  zenith). 
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Figure  2 


Let  Z  denote  the  zenith  distance  of  the  satellite,  and  y  the 
geocentric  angle  between  the  direction  to  zenith  and  direction 
to  the  SN  (Figure  1).  From  Figure  1  we  get  the  relation 

•  Jt 

«ta«  "  «teU-y)  (29  ) 

Solving  the  equation  with  respect  to  Z  we  get 

(30) 

•“y-— 

where : 

R  -  average  radius  of  Earth 
a  -  distance  of  SN  from  the  center  of  Earth 


For  an  SN  which  is  on  the  horizon  (Z  *  90°)  we  have 
according  to  (29) 


R 

eaayt  -  — 

(3D 

As  an  example,  assuming  R  *  6365  km  and  a  *  7^40  km  (the  value 
of  the  3emiJajor  axis  for  SN-12,  which  is  equal  approximately 
to  the  average  value  of  semimajor  axis  NSS  considered  in 
Tables  2  and  3)  we  get 

y* (32) 


If  the  SN  reaches  apogee  at  zenith,  then  the  maximum  stay  time 

At  of  a  given  SN  on  the  horizon  is  approximately 
m3  x 
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Jl„ 


»  — M 


(33) 


since  Au  *  2y 
Prom  the  data  for  SN-12  we  get  At 


h  and  we  have  agreement  with  (22) 


=  l8m45. 


max 


/215 


The  time  interval  which  elapses  between  the  rising 
(appearance)  of  a  given  SN  and  its  transit  through  zenith  is 
1/2  At. 


max 


The  above  values  of  vh  and  At  calculated  for  SN-12,  and 
the  conclusions,  are  also  valid  for  the  remaining  SN  of  the 
TRANSIT  system,  since  the  values  of  orbital  elements  as  such 
a,  e,  i,  n  are  very  close. 


The  time  interval  sutiable  for  Doppler  measurements  is 
shorter  than  the  maximum  staying  time  of  the  SN  above  the 
horizon,  since  In  computations  based  on  Integrated  Doppler 
observations  no  near-horizon  and  near-zenith  measurements  are 
used.  The  maximum  number  of  integrated  Doppler  type  measure¬ 
ments  -  taken  during  one  flight  of  an  SN  over  a  2-minute  time 
interval  does  not  exceed  8-9. 


In  turn,  we  shall  discuss  conditions  of  visibility  of  an 
SN  depending  on  the  geographic  latitude  of  the  observer. 
According  to  (16)  during  one  pass  of  an  SN,  its  orbit  will 
move  In  the  plane  of  equator  to  the  West.  For  example,  for 
SN-12  we  get 


Taking  this  into  consideration  and  the  value  (32)  we  see  that 
for  near-equatorial  regions,  the  SN  passes  over  the  horizon  of 
selected  point  only  during  2-3  consecutive  passes  -  once  on 
passing  through  the  equator  from  the  Southern  to  Northern 
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hemisphere  (ascending  orbit),  and  the  second  time  half  a 
day  later  on  passing  from  the  Northern  to  Southern  hemisphere 
(descending  orbit).  In  near-polar  regions  the  SN  will  be 
over  the  horizon  during  each  flight.  The  number  of  flights 
of  the  SN  per  day,  useful  for  geodetic  and  navigational  work, 
increases  with  geographic  latitude  of  the  observer  from  4-6 
in  areas  near  the  equator  to  12-14  in  near-polar  regions. 


For  an  observer  at  a  given  latitude  $0  the  SN  will  be 
over  the  horizon  only  when  the  longitude  of  the  ascending  node 

3.  h 

A  is  close  to  A  A  -  see  equations  (27)  and  (28).  The 
w  w  w 

difference  should  not  be  greater  than  the  value  of  &A  from 

m 

equation 


•in  J  A.,  = 


co* 


(35) 


This  equation  follows  from  the  spherical  triangle  OPS  (Figure  2). 
In  the  limiting  case  for  a  given  SN  reaching  apogee  on  the 
horizon  we  have  y  *yh  ,  J  l -Aim  .  The  internal  angle  at  the 
apexes  is  a  right  angle. 


Table  5 


..  ! 

Mm 

«•  ! 

Jim 

r  ! 

I 

HI 

sr 

M*7 

lr  i 

lit 

**•  ! 

M'S 

sr 

sr  i 

77*7 

lo- 

3T7 

srs 

•n 

40-  * 

4TS 

srts 

frt 

Table  5  gives  values  of  3Xm  calculated  according  to  (35) 
after  accepting  the  value  (32).  It  follows  also  from  equation 
(35)  that  for  AAm  *  rj 0°  we  have  Yh  *  90°  -  hence  for 
points  satisfying  the  condition 


the  SN  will  be  over  the  horizon  of  a  given  observer  during  each 
flight. 
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From  triangle  OPS  we  also  get  the  relation 


- 


tfy. 

cotd  Xk 


(37) 


where : 


=  (37a) 

iC  ic 

<J>g.  Xg  -  geographic  latitude  and  longitude  of  the 
subsatellite  point  at  the  time  of  apogee  of  a  given  SN. 


If  we  take  into  consideration  the  effect  of  the  rotational 
motion  of  the  Earth,  then  the  formula  for  geographic  latitude 
of  the  subsatellite  point  has  the  form 


i*/,  •> 


CU»  j  xk 


(6-3) 

|»-M 


(38) 


For  an  SN  culminating  rather  high  over  the  horizon 
(cos  a 1)  the  following  approximate  realtion  applies,  in 
conformity  with  (37)  or  (38) 


(39) 

SN 

For  aATRANSIT  satellite  at  the  time  of  apogee  the  geographic 
latitude  of  the  subsatellite  point  is  approximately  equal  to 
the  latitude  of  the  observer. 

For  example,  for  $0  *  52°  and  AX^  *  20°  the  SN  reaches 

apogee  at  the  zenith  distance  z  3  60°,  at  $  *  53°71  accord- 

s 

ing  to  (37)  or  at  <h  *  5h 0 2U  according  to  (38).  The  difference 

s  S 

<*>_  —  <t>o  “  2°24  is  covered  by  the  SN  during  4o  .  This  difference 
will  be  smaller  for  SN  reaching  apogee  higher  above  the  horizon 
and  larger  at  lower  altitudes.  In  computations  of  the  predict¬ 
ed  ephemeris  of  a  polar  SN  we  assume  that<pg  ■  4>0  -  see  equa¬ 
tions  (23),  (24)  and  (27). 


As  an  example,  let  us  consider  the  coordinates  of  a  point 
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in  the  area  of  Poland(A  =  20°,  4>  =  52°).  For  SN-12,  according 

to  (35)  we  obtain  AA  =  57°24.  Hence,  from  (27)  and  (28)  we 

m 

get 

A*-3rB6;  A‘--1S0OS® 

The  time  intervals  between  the  crossing  of  the  equator  by  a 
given  SN  and  apogee  are,  in  conformity  with  (24) 

IJ-I.-1S-4;  -  37-9  (Hi) 

Knowing  the  moment  of  passage  of  the  SN  through  the  node  we 
calculate  the  time  of  apogee,  on  the  basis  of  (4i).  The  SN 
will  fly  over  the  horizon  of  point  P  if  the  following  conditions 
are  fulfilled  for  the  longitude  of  the  ascending  node: 

a)  -33°4<aw<m°i  SN  moves  in  Soutn-North  direction 

(42) 

b)  ls'rss  «  A, <  26«aT£  SN  moves  in  North-South  direction 


Taking  the  values  of  (34)  and  (42)  we  see  that  a  given  SN  flies 
over  the  horizon  of  point  P  4-5  times  a  day  in  an  ascending 
orbit  and  similarly  on  a  descending  orbit.  Knowing  the  apogee 
point  tk  and  Aw,  we  calculate  from  (18)  the  longitude  of  the 
subsatellite  point  at  the  time  of  apogee  AK  and  (see  Figure  2) 


•in ft  -  cm 


(43) 


Knowing  the  geocentric  angle  r*  ,  we  calculate  from  equation 
(30)  the  zenith  distance  of  the  SN  at  the  time  of  apogee.  At 
the  ascending  or  descending  time  of  the  SN  we  have  r*r*. 
hence  from  the  cosine  formula  for  spherical  triangle  OS^S^ 
(Figure  2)  we  have 

eMS,S|  - 

«"ft  (44) 

I  0  I 

Since  S,  Sk  ■  (n-|w|At,  the  time  interval  which  elapses  between 
the  ascension  and  apogee,  or  apogee  and  descent  of  the  SN  is 
expressed  by  the  equation 

ti,  -  -J-> 

•-I*  (45) 
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Equations  (44)  and  (45)  are  approximate,  since  they  were 
deduced  on  the  assumption  that  AXk  =  AXEast  =  AXWest,  the 
angle  AX  does  not  change  during  the  motion  of  the  SN  over  the 
horizon  and  is  equal  to  its  value  at  apogee,  or  the  average 
value  during  the  flight  of  the  SN  over  the  horizon. 

On  the  basis  of  these  formulas  we  can  calculate  the  ephemeris  of  an  SN 
for  any  pcirt  on  the  surface  of  the  Earth  during  the  day.  ' 

On  the  next  day,  the  conditions  of  flight 
over  the  horizon  of  a  given  point  repeat  themselves  approxi- 
mately  after  12-13  passes  around  the  Earth.  Hence,  we  add  the 
values  of4AAw,  12  T  or  13  AXw,  2-3  T  to  the  corresponding  values 
of  ascending  Xw,  t  at  which  SN  appears  on  the  horizon. 

The  discussed  equations  can  be  used  to  prepare  computed 
ephemerides  for  SN  on  a  computer.  In  [3]  there  is  a  published 
program  in  Fortran  language  to  calculate  the  ephemerides  for 
an  SN  at  medium  latitudes  and  in  areas  near  the  equator  omitting 
the  polar  regions. 

SN 

The  data  necessary  to  predict  the  ephemerides  of^ TRANSIT  /2l6 
satellites  are  also  provided  by  the  Defense  Mapping  Agency, 
Hydrographic  Center,  Washington,  D.C.  20390  [3].  These  give 
the  values  of  Xw,  tw,  AX^,  T.  Table  6  gives  the  same  parameters 
calculated  for  6  navigation  satellites  for  December  17,  1975. 

They  were  calculated  on  the  basis  of  orbital  elements  given  in 
Table  3  and  for  SN-18  In  Table  2.  The  parameters  and  orbital 
elements  given  in  Tables  2  and  3  were  taken  from  Doppler 
measurements  of  a  TRANSIT  SN  made  by  the  Higher  Naval  Academy 
in  Gdynia. 

Table  6 


N.  SN 

1  «»  JTU)  * ! 

1  T  ! 

— isriw* 

io*»ss;» 

— 26*7M 

44*4704  1 

— lorssu 

147*0449 

—-24*017 

IM-SM4 

—  144*161$ 

104*7990 

—24*770 

449*069f 

■*. 78*24*3 

104*9764  | 

— 24*017 

J«*4*?620 

•  «**SWT 

107*4249 

— 26*029 

Jtl 

*  7<l434iJ 

105*619$ 

—24*474 
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5.  Graphic  Method  of  Predicting  SN  Flights 

Predictions  for  the  flight  conditions  of  SN  over  a  given 
point  on  the  surface  of  Earth  can  be  considerably  accelerated 
by  the  use  of  a  graphic  method.  Figure  3  shows  a  contour 
map  in  Mercator  projection  with  the  point  P(  20°,  $  =  52°  - 

central  Poland).  At  this  point  a  net  is  drawn  for  equal 
altitudes  (almucantars )  for  the  following  horizontal  elevations 
0°(horizon),  7°5,  30°  and  60°,  and  of  azimuths  (verticals) 
marked  in  intervals  of  30°  in  the  range  0°-360°.  This  net  was 
calculated  for  a  SN  circling  at  a  distance  of  a  =  7440  km  from 
the  center  of  Earth.  On  a  transparent  paper  and  in  the  same 
projection,  the  path  of  several  consecutive  SN  passes  (for 
ascending  and  descending  orbits)  is  drawn.  In  this  case  it  is 
sufficient  to  know  only  the  values  of  Xw,  tw,  T  from  one  (any) 
flight  on  a  day  to  be  able  to  determine  visibility  conditions 
of  a  SN  at  any  point  on  the  surface  of  Earth.  As  an  illustration, 
4  consecutive  flights  of  SN-12  (in  ascending  and  descending 
orbits)  visible  over  the  horizon  in  Poland  on  17  December  1975 
were  drawn  in  Figure  3.  The  graphs  were  made  on  the  basis 
of  data  given  in  Tables  3  and  6.  Assigning  No.  1  to  an  orbit 
whose  ascension  is  defined  by  the  parameters  t  =  263m6673(TU) , 

Xw  -  -128°109,  the  following  passes  will  be  visible  at  point 
P:  0,  1,  2,  3  (descending  orbit)  and  7,  8,  9,  10  (ascending 
branch  of  the  orbit,  marked  in  Figure  3  by  dots  every  2m) . 

In  order  to  rapidly  determine  the  moment  of  time  at  which 

the  SN  will  be  over  a  defined  point  on  the  surface  of  the 

Earth  (for  determination  of  the  time  of  apogee  t.  ) ,  the  position 

in  ** 

of  the  SN  was  marked  every  2  ,  counting  from  zero  at  the 
ascending  node.  Thus,  for  example,  for  pass  No.  2  over  the 
horizon  at  point  P  we  get  t*  *  tw  +  T  +  38m =  6h48m2.  For  pass 
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Figure  3 
(Rownik 


L 


No.  8  we  obtain  =  t  +  7T  +  l6m  =17h5m4.  During  passes 
No.s  0,  1  the  satellite  moves  on  the  Eastern  side  of  the  sky, 
and  during  passes  No.s  2,  3  -  on  the  Western  side.  In  these 
A  cases  the  ascension  of  the  SN  is  on  the  Northern  side  of 
the  horizon,  and  the  descent  -  on  the  Southern  side.  In  turn, 
during  passes  No.s  7,  8,  9,  10  the  ascension  is  on  the  Southern 
side,  and  the  descent  -  on  the  Northern  side  of  the  horizon. 

The  passes  No.s  7,  8  take  place  on  the  Eastern  side,  and 
No.s  9,  10  -  on  the  Western  side  of  the  horizon. 
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